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Abstract 

Fear behaviors and fear memories in rodents have been traditionally assessed by the amount of freezing upon the 
presentation of conditioned cues or unconditioned stimuli. However, many experiences, such as encountering earthquakes 
or accidental fall from tree branches, may produce long-lasting fear memories but are behaviorally difficult to measure using 
freezing parameters. Here, we have examined changes in heartbeat interval dynamics as physiological readout for assessing 
fearful reactions as mice were subjected to sudden air puff, free-fall drop inside a small elevator, and a laboratory-version 
earthquake. We showed that these fearful events rapidly increased heart rate (HR) with simultaneous reduction of heart rate 
variability (HRV). Cardiac changes can be further analyzed in details by measuring three distinct phases: namely, the rapid 
rising phase in HR, the maximum plateau phase during which HRV is greatly decreased, and the recovery phase during 
which HR gradually recovers to baseline values. We showed that durations of the maximum plateau phase and HR recovery 
speed were quite sensitive to habituation over repeated trials. Moreover, we have developed the fear resistance index based 
on specific cardiac response features. We demonstrated that the fear resistance index remained largely consistent across 
distinct fearful events in a given animal, thereby enabling us to compare and rank individual mouse's fear responsiveness 
among the group. Therefore, the fear resistance index described here can represent a useful parameter for measuring 
personality traits or individual differences in stress-susceptibility in both wild-type mice and post-traumatic stress disorder 
(PTSD) models. 
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Introduction 

Emotionally charged episodic events can often produce robust 
and long-lasting memories even upon a single exposure [1-3]. 
Extensive behavioral, lesion and brain imaging studies have shown 
that emotion can facilitate the encoding, consolidation and 
retrieval of memory [4-7]. In addition, emotional stimuli, 
especially the fearful/aversive stimuli, typically evoke innate 
unconditioned responses [8-11]. Therefore, many fearful stimuli 
have been used as unconditioned stimuli in classical fear 
conditioning to study the neural mechanism of learning and 
memory [12-17] and are useful probes for investigation of 
emotions [18,19]. 

In rodents, freezing behavior has been widely used to assess 
emotional responses elicited by fearful stimuli. Freezing is defined 
as the absence of body movement except for breathing in rodents 
[9, 1 0] . The other type of response, startle response, which consists 
of a contraction of the skeletal, facial and neck muscles, is also used 
to assess fear [20-23]. In humans, the eye blink behavior as a 
startie response to air pufF to eyes is typically measured by 



electromyography (EMG) of the orbicularis oculi muscle [24] . In 
rats and mice, stabUimeter and accelerometer devices are also used 
to record the amplitude of the starde response to fearful events 
[25-27]. 

Recently, we have used our large-scale neural ensemble 
recording technique to monitor the real-time encoding patterns 
of fearful experiences in hippocampal CAl region and anterior 
cingulate cortex (ACC) of freely behaving mice [28—30]. We have 
shown that dilTerent fearful events, such as air puff, free fall, and 
earthquake-like shake, evoked diverse responses in CAl and ACC 
neurons, consequently permitting quantitative measurement of the 
formation of real-time neural ensemble traces of these fearful 
experiences in the mouse brain [28-30]. In addition, we have 
shown that earthquake and free fall can evoke robust suppression- 
and-offset-rebound excitation of the ventral tegmental area (VTA) 
putative dopamine neurons in freely behaving mice [31], 
indicating the termination of aversive or fearful events are strong 
drive for increasinging dopamine activity in the brain. Unfortu- 
nately, the traditional behavioral responses, such as the freezing 
behavior, can't be used to quantitatively measure fearfulness 
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evoked by earthquake, free fall, or air puff because animals often 
move and exhibit various escape behaviors during or immediately 
after these events. Similarly, the transient starde responses [2,26], 
measured by stabilimeter and accelerometer devices, are not 
suitable to provide continuous assessments under such experimen- 
tal paradigms. 

These limitations have prompted us to examine autonomic 
responses, such as changes of heartbeat, as a method reported in 
the literature to assess behavioral changes in rats [32-36]. HR and 
HRV can be determined from electrocardiogram (ECG) gener- 
ated by the sinoatrial node. It has been shown that changes of HR 
and HRV are correlated with emotional changes in animals, such 
as fear in both mice and rats [37-39]. In addition, several early 
works in rats also reported HR changes to startUng acoustic sound 
or air puff stimuli [33-35,40,41]. Consistent with these findings, 
our recent mouse study has further demonstrated that fear 
conditioning can trigger strong changes in HR and HRV in 
freely behaving mice [42]. In particular, we have further shown 
that HRV is highly correlated with the performances in both 1-h 
short-term and 24-h long-term fear memory tests [42] . 

In the present study, we set out to examine whether and how 
emotionally fearful events, such as earthquake, free fall, and air 
puff, would trigger cardiac changes in freely bcha\'ing mice. We 
have systematically and quantitatively characterized the various 
stages of HR responses to air puff, free fall, and earthquake-like 
cage-shake in freely behaving mice. Moreover, we have examined 
whether and how multiple trials or exposures to such fearful events 
produce habituation-like changes in HR and HRV dynamics. 
Finally, we have employed the fear resistance index to describe 
and compare the summed autonomic responses across individual 
mice. 

Results 

Three Distinct Heart Rate Response Phases to Fearful 
Stimuli in Freely Behaving Mice 

To investigate whether and how air puff, free fall, and 
earthquake alter the emotional states, we set out to measure the 
changes in HR and HRV in freely behaving mice. We implanted a 
pair of electrodes subcutaneously into the chest of these mice for 
ECG recording using the procedures previously described [42]. 
Data from seven mice, in which the ECG signals were stably 
recorded, were used for the current analyses. The ECG data 
showed that air puff, free fall, and earthquake stimuh all triggered 
significant increases in HR (Figure lA & B). 

Moreover, we found that there were three distinct phases of HR 
responses to these stimuli: the rising phase, the maximum plateau 
phase and the recovery phase (Figure IC). The rising phase was 
between the onset of stimulus and the start of HR plateau; the 
plateau phase was the stable period of maximum HR during 
which HRV was greatly reduced; and the recovery phase was the 
duration for HR gradually coming back to the basal level (see 
methods for details). 

Changes in Heart Rate and Heart Rate Variability in 
Response to Sudden Air Puff 

To study how the dynamics of HR and HRV would respond to 
air puff, we delivered a sudden air blow (10 psi for 0.5 s) to the 
back of the mouse, and repeated for seven timc-s with 2-5 min 
randomized intervals between each air puff. In previous startle (air 
puff) studies in rats, both the tachycardia and bradycardia have 
been reported [33,40,41]. However, we only observed the 
tachycardia, a robust increase of instant HR after the onset of 
air puff stimulus in mice (Figure 2A, 1*' Trial). It is noted that 



repeated air puff to the animals reduc(-d its effectiveness in eliciting 
cardiac changes, indicating the habituation effect (Figure 2A, 4''' 
and 7* Trials). In order to quantitatively measure the real-time 
HRV, we performed a nonlinear dynamic technique, Poincare 
plot analysis, which can be applied for visualization of the R-R 
interval fluctuation and assessment of the dynamics of HRV. 
Poincare plot of R-R intervals revealed prominent increase in 
heartbeat regularity (or reduction in HRV) during the 30-s post- 
stimulus period (Figure 2B, the left panel indicates the response to 
Trial). The Poincare plot analysis revealed the habituation 
effect over seven trials (Figure 2B, the middle and right panels 
showed cardiac responses to the 4''' and 7''' Trials, respectively). 
The mean HR responses from the seven mice showed that HR 
increased significantiy from the basal level of 494±27 bpm to 
690±24 bpm to the first air puff (P<0.001) (Figure 2C). More- 
over, analysis of the averaged coefficient of variation (C V) of R-R 
intervals, measurements of HRV, has revealed the noticeable 
decrease after the 0.5-s air puff (Figure 2D). In addition, the root 
mean square of the successive differences (RMSSD) of R-R time 
intervals also confirmed decrease in HRV (Figure 2E). 

To provide further quantitative measurement of the habituation 
effects on the dynamic HR responses to air puff, we evaluated five 
distinct components of cardiac responses: namely, (1) the rising 
speed (of the rising phase) which characterizes the speed of the 
elevation of HR upon stimulation; (2) time duration of the 
maximum plateau; (3) HR of the maximum plateau; (4) C V of the 
maximum plateau; and (5) reccn ery speed (of the recover)' phase). 
Interestingly, our analysis suggests that the rising speed did not 
change significantly throughout seven trials (Figure 3A, 
= 0.14, /*>(). 05). The maximum plateau duration time showed a 
significant decrease from the second trial and there afterwards in 
comparison to that of the first trial (Figure 3B, i^g 33 = 6.55, P< 
0.05, P<0.01, P<0.001). Yet, HR was consistenfly increased over 
the baseline after receiving air puff across all seven trials 
(Figure 3C, Fj^ 42 = 14.08, P<0.001), while CV in the maximum 
plateau phase exhibited dramatic and consistent decrease over the 
baseline after receiving air puff (Figure 3D, i^y 42 = 5.50, P<0.01). 
By contrast, the recovery speed showed a statistically significant 
increasing trend as air puff was repeated (Figure 3E, 35 = 4.71, 
P<0.05, P<0.01). These observations showed that distinct cardiac 
parameters exhibited different responses to repeated air puff 
stimulation. In particular, the time duration of the maximum 
plateau phase and the HR recovery speed were quite sensitive for 
detecting habituation effects over multiple trials. 

Changes in Heart Rate and Heart Rate Variability in 
Response to Free Fall 

A sudden free fall can be an intense fearful event which mimics 
the tower-of-terror experience in Disney world for humans or the 
accidental fall from a tree branch for small animals. Our recent 
in vivo neural recording in mice revealed that free fall triggered 
dynamic firing changes in many neurons in the CAl region of the 
mouse hippocampus [28] and also triggered firing changes in three 
types of the VTA dopaminergic neurons [31]. However, changes 
in HR dynamics during free fall in rodents have not been 
investigated. Thus, to measure changes in cardiac responses to free 
fall (or elevator drop), we placed the mouse in a small chamber 
which was 1 3 cm off the ground. In order to minimize the effect of 
novel environment on basal HR [42] , mice were habituated to this 
chamber for three days (1 h each day). On the day of the 
experiments, we delivered a computerized signal to initiate the 
drop, and both HR and HRV were measured. The trials were 
repeated for seven times with 2-5 min randomized time intervals 
between each drop trial. We found that the instant HR increased 
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Figure 1. ECG recording and tKiree distinct pKiases of HR responses to fearful stimuli. (A) Illustration of three fearful stimuli paradigms: 
namely, air puff, a sudden air blow (10 psi, 0.5 s) to the back of the animal; free fall, an abrupt drop of the animal inside a small chamber (13 cm high); and 
earthquake, an unexpected earthquake-like shaking to the animal (3000 rpm, 0.5 s). (B) Examples of ECG recording during air puff free fall and earthquake 
stimuli. The red dots indicate the peaks of the R-wave. The blue triangles indicate the onsets of each stimulus. Please note the significant shortened R-R 
intervals after the onset of all the three types of stimuli. Although noise was introduced during the presence of the stimuli, the custom-written Matlab 
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program can still precisely extract the timestamps of the R-wave peaks. Scales: 0.5 mV. (C) Schematic representation of the three distinct phases of HR 
responses to fearful stimuli: first (t1, the purple triangle), the rising phase; second (t2, the red rectangle), the maximum plateau phase; and third (t3, the 
cyan triangle), the recovery phase. 
doi:1 0.1 371 /journal.pone.0093270.g001 



in response to drop episodes (Figure 4A). Typically, the cardiac 
responses became smaller over repeated trials (Figure 4A). The 
Poincare plot analysis also showed dynamic patterns in HRV 



decrease immediately and remained so for a certain period of time 
(Figure 4B). As shown in the Figure 4B, the initial 30 s (in red) in 
the first trial showed the tight R-R time interval distribution in this 



Air puff 




Time (s) 

Figure 2. Response of HR dynamics to air puff. (A) Instant HR responses of an individual mouse during three trials of air puff stimuli, trial #1, 
trial #4 and trial #7. The red vertical dashed lines indicate the onset of the air puff stimuli (10 psi, 0.5 s). Robust increase of instant HR after stimulus 
onset was noted at trial #1. (B) Poincare plots of the same mouse's R-R intervals of 30-s pre-stimulus (grey line) and 30-s post-stimulus (red line) 
period in trial #1, trial #4 and trial #7. Prominent increase in heartbeat regularity during the 30-s post-stimulus period was noted in trial #1. (C) The 
average HR response to air puff (l" trial). (D) The average CV of R-R intervals upon l" trial of air puff. (E) The average RMSSD of R-R intervals upon l" 
trial of air puff. The average curves of 30-s pre-stimulus period were plotted as grey lines, whereas the average curves of 90-s post-stimulus period 
were plotted as red lines. The shadows around the average curves indicate error bars (SEM); n = 7. 
doi:1 0.1 371/journal.pone.0093270.g002 
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Figure 3. Effects of habituation on HR dynamics upon multiple trials of air puff. (A) The average rising speed showed no significant 
difference during seven trials of air puff stimuli, n = 7. (B) The plateau durations of seven air puff stimuli exhibited a significant decreasing trend from 
the second trial and there afterwards. n = 7; *P<0.05, **P<0.01, ***P<0.001, one-way repeated measures ANOVA and Dunnett's multiple 
comparisons test. (C-D) In comparison to the pre-stimulus state, the significant changes of HR and CV during the maximum plateau maintained 
throughout seven trials. n = 7; **P<0.01, ***P<0.001, one-way repeated measures ANOVA, Tukey post hoc test and Dunnett's multiple comparisons 
test. (E) The recovery speed across seven air puff stimuli showed a significant increase as air puff was repeated, n = 7; *P<0.05, **P<0.01 , one-way 
repeated measures ANOVA and Dunnett's multiple comparisons test. All data are plotted as mean ± SEM (error bars). 
doi:1 0.1 371 /journal.pone.0093270.g003 



representative mouse (the left panel). However, the same 30-s 
periods from the 4* and 7* drop showed more chfiFused 
distribution (the middle and right panels, respectively). As a 
group, HR showed significant increases upon free fall drop. For 
example, the averaged HR was up-regulated from 567 ±28 bpm 
at the basal level to 731 ±9 bpm, (P<0.001) upon the T' drop, and 
then gradually decreased over the next 90 s (Figure 4C). For the 
same period of time, the averaged coefficient of variation (CV) of 
R-R intervals showed a decrease upon the drop (Figure 4D). 
Similarly, the root mean square of the successive differences of R- 
R time intervals also showed a decrease in HRV (Figure 4E). 

To examine in details the effects of multiple free fall events on 
HR changes, we evaluated the rising speed, time durations of the 
maximum plateau, HR and CV of the maximum plateau, and 
recovery speed over multiple drop trials. Although there was no 
habituation trend in HR acceleration (Figure 5A, F(j 155 = 1.19, P> 
0.05), the time durations of the maximum plateau phase showed 
significant decreases over seven trials (Figure 5B, Ff, 35 = 3.72, P< 



0.05, P<0.01, P<0.001). The elevation of HR and reduction of 
CV continued throughout seven trials (Figure 5C & D, 42 
= 24.60, P<0.001, i?7, 42=15.44, P<0.001). In addition, as a 
group, the recovery speed of the HR back to the basal level also 
significandy increased by the 5* trial and afterwards (Figure 5E, 
i^s, 36 = 5.43, P<0.05, P<0.001). Therefore, these detailed 
analyses revealed that while the HR and HRV continually showed 
significant differences over the basal levels, the time durations of 
the maximum plateau and the recovery speed were changing as 
drop was repeated over seven trials. 

Changes in Heart Rate and Heart Rate Variability in 
Response to Earthquake 

To examine how earthquake events alter the dynamic HR 
responses, we placed mice in a small circular chamber which was 
fixed on top of a vortex machine, and earthquake event was 
triggered by a computerized program. Again, the mice underwent 
1-hr habituation session each day for three days to minimize 
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Figure 4. Response of HR dynamics to free fall. (A) Instant MR responses of an individual mouse during three free fall stimuli, trial #1, trial #4 
and trial #7. The red vertical dashed lines indicate the onset of the free fall stimuli (13 cm high). (B) Poincare plots of the same mouse's R-R intervals 
of 30-s pre-stimulus (grey line) and 30-s post-stimulus (red line) period in trial #1, trial #4 and trial #7. Immediate decrease in dynamic patterns in 
HRV after stimulus onset was noted at trial #1 , but more diffused distribution in the 30-s post-stimulus period in trial #4 and trial #7. (C) The average 
HR response to free fall (l" trial). (D) The average CV of R-R intervals upon l" trial of free fall. (E) The average RMSSD of R-R intervals upon l'^' trial of 
free fall. The average curves of 30-s pre-stimulus period were plotted as grey lines, whereas the average curves of 90-s post-stimulus period were 
plotted as red lines. The shadows around the average curves indicate error bars (SEiVl); n = 7. 
doi:10.1371/journal.pone.0093270.g004 



effects of the novel chamber. On the day of experiments, the mice 
were subjected to seven trials of earthquakes with 2-5 min of 
randomized time intervals between each trial. The instant HR 
showed the rapid tachycardia to earthquake. Responses of a 
representative mouse to the P', 4* and 7* shake trials are shown 
in Figure 6 A (from the left to right panels). Although the 
earthquake occurred for only 500 ms, the duration of the 



maximum plateau phase could last up to 30 s (i.e. the first 
earthquake trial). The reduction in HRV in response to 
earthquake was further evident from Poincare plot analysis of 
the given individual mouse (Figure 6B, the initial 30-s R-R time 
distribution plot from the f' Trial, red fine). On average, HR 
increased from basal level of 51 8± 18 bpm to 725±15 bpm upon 
the 1'^ earthquake (P<0.001) (Figure 6C). Moreover, the averaged 
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Figure 5. Effects of habituation on HR dynamics upon multiple trials of free fall. (A) The average rising speed showed no significant 
difference during seven trials of free fall. n = 7. (B) The time duration of the maximum plateau displayed a significant decrease over trials. n = 7; *P< 
0.05, **P<0.01, one-way repeated measures ANOVA and Dunnett's multiple comparisons test. (C-D) In comparison to the pre-stimulus state, the 
significant changes of HR and CV during the maximum plateau maintained throughout seven trials. n = 7; ***P<0.001, one-way repeated measures 
ANOVA, Tukey post hoc test and Dunnett's multiple comparisons test. (E) The recovery speed across seven free fall stimuli showed a significant 
tendency in becoming faster by the 5* trial and afterwards, n = 7; *P<0.05, ***P<0.001, one-way repeated measures ANOVA and Dunnett's multiple 
comparisons test. All data are plotted as mean ± SEM (error bars). 
doi:1 0.1 371 /journal.pone.0093270.g005 



CV and RMSSD of R-R intervals also showed reductions in HRV 
(Figure 6D & E). 

Furthermore, we have systematically examined the various 
parameters on dynamic HR responses upon multiple trials of 
earthquake. The slope of HR increase maintained similar across 
trials (Figure 7A, F^^ 35 = 0.79, P>0.05). However, time durations 
for the maximum plateau phase decreased significantly from the 
second earthquake trial and then afterwards (Figure 7B, F^^ 
= 7.71, P<0.05, P<0.01, P<0.001). Interestingly, the HR and 
CV of maximum plateau still existed significant responses 
throughout seven trials (Figure 70 & D, F^^ 42 = 66.70, P<0.001, 
Fj^ 42 = 20.46, P<0.001). On the other hand, the recovery speed 
exhibited a tendency in becoming quicker over trials, with 
statistically difiFerences at the 5* trial and afterwards (Figure 7E, 
i?6, .■56= 3.59, P<0.05, P<0.01.). 



Measurement of Summed Fear Responses to Different 
Events among Individual Animals 

Under many circumstances, it may be highly desirable to 
compare the fear responses among different animals. Moreover, 
we are also interested in asking whether an individual animal's fear 
response to a given fearful event can be used to predict the 
animal's responses to other fearful events. In order to address these 
questions, we set out to develop and calculate the fear resistance 
index by combining five individual HR parameters: (1) the rising 
speed (of the rising phase); (2) time duration of the maximum 
plateau; (3) HR of the maximum plateau; (4) CV of the maximum 
plateau; and (5) recovery speed (of the recovery phase). The fear 
resistance index, calculated from comprehensive evaluation index 
by using the above five factors based on the Shannon's entropy 
method [43-45], is the quantitative measurement of summed 
responses to fearful events. The higher the fear resistance index, 
the lower the HR response magnitude of mice to fearful challenges 
(see methods for details). 
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Figure 6. Response of HR dynamics to earthqual<e. (A) Instant MR responses of an individual mouse during three trials of earthquake stimuli, 
trial #1, trial #4 and trial #7. The red vertical dashed lines indicate the onset of the earthquake (0.5 s). (B) Poincare plots of the same mouse's R-R 
intervals of 30-s pre-stimulus (grey line) and 30-s post-stimulus (red line) period in trial #1, trial #4 and trial #7. Note the reduction of HRV in 
response to earthquake at trial #1, and such reduction effect became weaker at trial #4 and trial #7. (C) The average HR response to earthquake (1^' 
trial). (D) The average CV of R-R intervals upon l" trial of earthquake. (E) The average RMSSD of R-R intervals upon l" trial of earthquake. The average 
curves of 30-s pre-stimulus period were plotted as grey lines, whereas the average curves of 90-s post-stimulus period were plotted as red lines. The 
shadows around the average curves indicate error bars (SEM); n = 7. 
doi:1 0.1 371 /journal.pone.0093270.g006 



We first compared the fear resistance indexes during air puff, 
free faU and earthquake with those during neutral tone. As 
expected, the fear resistance indexes during fearful events were 
significantly lower than those during neutral tone (Figure 8, i^J_ 23 
= 10.44, P<0.01), therefore providing evidence for the validity of 
such index. 



Next, we calculated the fear resistance indexes of the mice 
during air puff. Our analysis showed that Mouse #1, #2, and 
#7 were ranked as the top three (in order) for exhibiting the 
strongest HR responses to air puff, whereas Mouse #6 had the 
least response to sudden air puff and Mouse #5 the second least 
response (Figure 9A, the left panel). At the same time, Mouse #3 
and #4 belonged to the middle of the group (Figure 9A, the left 
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Figure 7. Effects of habituation on HR dynamics upon multiple trials of earthquake. (A) The average rising speed showed no significant 
difference during seven trials of earthqual<e stimuli. n = 7. (B) The plateau durations of seven earthquake stimuli decreased significantly from the 
second trial and afterwards, n = 7; *P<0.05, **P<0.01, ***P<0.001, one-way repeated measures ANOVA and Dunnett's multiple comparisons test. (C- 
D) In comparison to the pre-stimulus state, the significant changes of HR and CV during the maximum plateau maintained throughout seven trials, 
n = 7; ***P<0.001, one-way repeated measures ANOVA, Tukey posf hoc test and Dunnett's multiple comparisons test. (E) The recovery speed across 
seven earthquake stimuli exhibited a significant increasing trend over trials. n = 7; *P<0.05, **P<0.01, one-way repeated measures ANOVA and 
Dunnett's multiple comparisons test. All data are plotted as mean ± SEM (error bars). 
doi:1 0.1 371 /journal.pone.0093270.g007 
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Figure 8. Comparison of fear resistance indexes between three 
types of fearful events and neutral tone. The fear resistance 
indexes during fearful events (n = 7) were significantly lower than those 
during neutral tone (n =6). **P<0.01, one-way ANOVA and Tukey posf 
hoc test. All data are plotted as mean ± SEM (error bars). 
doi:1 0.1 371/journal.pone.0093270.g008 



panel, resistance indexes as follows: mouse #1< #2< #7< #3< 
#4< #5< #6). Second, we calculated the summed responses to 
free fall. Interestingly, the fear resistance indexes obtained from 
free-fall drop experiments also showed the similar rank order with 
Mouse #1, #7, and #2 as the top three mice for exhibiting the 
strongest HR-related dynamics, and whereas mouse #4 and 
^3 in the middle and Mouse #6 and ^5 had the least fearful 
responses (Figure 9A, the middle panel, mouse #1< #7< #2< 
#4< #3< #6< #5). Finally, we further calculated the summed 
responses to earthquake. We also found that the same group of 
mice exhibited the identical rank order of their fear resistance 
indexes for earthquake in comparison to those of air puff 
(Figure 9A, the right panel, mouse #1< #2< #7< #3< #4< 
#5< #6). These results suggest that fear responses exhibited by 
each individual mouse tended to be rather consistent when 
exposed to different fearful events, thereby indicating the possible 
epigenetic factors in defining and controlling fear responses in a 
given individual. 
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Figure 9. Consistency of cardiac responses to three different types of fearful stimuli in individual mice. (A) For each type of fearful 
stimuli (air puff, free fall, and earthquake, respectively), the fear resistance index (FRI) was developed and calculated from comprehensive evaluation 
index by combining five individual HR parameters. Individual mice were ranked in the order of the fear resistance index (air puff, the left panel; free 
fall, the middle panel; earthquake, the right panel). The fear resistance index of each mouse revealed the high consistency of cardiac responses to the 
three types of fearful stimuli in individual mice. (B) The regression analysis of fear resistance index (FRI) in seven mice between air puff and free fall 
(the left panel), between free fall and earthquake (the middle panel), and between earthquake and air puff (the right panel) was carried out 
respectively. Significant correlations were observed across all the three different types of fearful stimuli. 
doi:1 0.1 371 /journal.pone.0093270.g009 



Further assessing this consistency of animals' summed fearful- 
ness across air puff, free fall, and earthquake, we found the 
summed fear responses from these seven mice during air puff and 
free fall exhibited significant correlation (r = 0.865, P = 0.0 11 9) 
(Figure 9B, the left panel). Similarly, the correlation between the 
fear resistance indexes for free fall responses and earthquake 
responses was also statistically significant (Figure 9B, the middle 
panel, r = 0.838, P = 0.0 185). In addition, the regression analysis 
also showed the tight correlation between fear resistance indexes to 
earthquake and air puff (Figure 9B, the right panel, r = 0.978, 
/" = 0.0001). These analyses again suggest the consistent cardiac 
responses of individual mice to different fearful stimuli. 

Discussion 

Our above experiments provided a detailed characterization of 
changes in HR and HRV, and recovery speed in freely behaving 
mice subjected to three different types of fearful events. Our use of 
free fall paradigm has been prompted by our desire to understand 
how and why Bungee jump or the Tower of Terror ride at Disney 
World would create long-lasting memory in humans. While such 
experiences can be highly fearful, our recent study in VTA 
revealed that some dopamine neurons exhibited suppression- 
rebound excitation dynamics. This may explain why such events 
can be paradoxically rewarding [31]. The stimulation paradigms, 
such as free fall and earthquake, which are quite relevant to both 



humans and animals, have also expanded experimental protocols 
for the study of fear memory processing that was often traditionally 
relied on fear conditioning paradigms or air puff. Over the past 
several years, we have investigated in vivo neural activity in several 
mouse brain regions including the CAl, ACC and VTA during 
free fall, earthquake, and air puff [28-31]. We have shown that 
our novel behavioral protocols are quite useful to enrich our 
understanding how the hippocampus and other brain regions 
encode fear memories [29,46]. Three major findings emerged 
from our present study: First, earthquake, free fall, and air puff all 
evoked significant changes in cardiac responses such as increased 
HR and decreased HRV. It is stated that HRV could be regarded 
as a reflection of moment-to-moment interplay between the two 
countervailing branches of the autonomic nervous system, namely 
an excitatory sympathetic component and an inhibitory parasym- 
pathetic component [47,48]. Thus changes in HRV were 
associated with changes in the regulatory output via efferent 
nerve traffic of the autonomic nervous system to the sinoatrial 
node [38,49]. In the literature, both increased HR and decreased 
HRV were typically observed in mice during acquisition as well as 
contextual and auditory cued retention phase of the aversive 
conditioning, a form of associative learning [37,38,42]. Such effect 
on HR during tone recall phase could be inhibited by the 
nonspecific beta-adrenergic antagonist sotalol, indicating a strong 
sympathetic activation [37] . Moreover, several lines of evidence 
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suggested that reduced HRV but not HR could serve as a more 
sensitive marker of enhanced fear response [42,47] . 

While we are not aware of any study conducted to measure HR 
and HRV during free fall and earthquake-like shake, several 
previous studies have reported the autonomic responses evoked by 
air puff and acoustic startle [32-36,40,41,50]. It has been reported 
that alerting stimuli such as acoustic startle or air puff stimuli could 
elicit tachycardia in freely moving mice [26,47], as well as in 
unrestrained rats [26]. However, some studies showed bradycardia 
in restrained rats when they were subjected to air puff, and such 
cardiac deceleration could be habituated over trials [33,40,41]. 
The cardiac responses to such alerting stimuli could reflect the 
dynamic changes in sympathetic and parasympathetic modulation 
under different behavioral states. In comparison, our finding in 
freely behaving mice showed a robust and consistent tachycardia 
to all the three fearful stimuli, but no bradycardia was observed. It 
is noteworthy to point out that upon the earthquake-like shake and 
drop, mice tended to move or escape. Thus, the traditional 
behavioral parameters such as the amount of freezing were not 
suitable for measuring fearful responses under these situations. 
The demonstration of changes in HR and HRV in mice should 
provide strong evidence for the notion that earthquake and free 
fcdl indeed represent fearful experiences to these rodents, just as 
these events to humans. This may, in turn, support the idea that 
these events can be and should be meaningful, not only 
ecologically, but also experimentally, to the study of fear memories 
beyond the traditional fear conditioning paradigms using context 
or cue associated with mild foot-shock. 

Second, to further understand how cardiac physiolog)' changes 
over repeated exposures to earthquake-like shake, free fall, or air 
puff to the back of the animals, we have divided cardiac dynamic 
responses into five distinct components, namely, the rising speed 
(of the rising phase), time duration of the maximum plateau, HR 
of the maximum plateau, CV of the maximum plateau, and 
recovery speed (of the recovery phase). Our systematic analyses 
have revealed that each of these five components contained useful 
information about the cardiac responses to these fearful events. 
These distinct components exhibited differential responses to 
habituation protocols. Specifically, the HR acceleration did not 
differ across trials. However, the HR and HRV of the maximum 
plateau phase showed significant differences upon the stimuli but 
remained relatively similar over all seven trials. Thus, these 
observations indicate that the slope of HR increase, increment of 
HR and decrement of HRV are important to reflect the cardiac 
response, but are not necessarily suitable to reflect the habituation 
effect. Furthermore, this may be different with the extinction of 
tachycardia triggered by conditioned stressor, during which the 
maximum HR decreased with repeated testing [5 1] . Interestingly, 
we found that time durations of the maximum plateau and the 
recovery speed were more sensitive to show such habituation 
effects. In addition, the highest baseline HR values and the largest 
baseline variation were noted during seven trials of free fall. The 
chamber needed to be raised vertically before each trial during 
free fall experiments, which could be the cause of the baseline 
fluctuation. Because the HR acceleration slowed down with higher 
HR [26], such fluctuation could lead to the variation in rising 
speed of HR (Figure 5A). Since the rising speed was one of the 
parameters for fear resistance index calculation, it should also 
result in the lower correlations in Figure 9 (left and middle panels). 

Third, we employed a Shannon's entropy method and 
developed a novel method to describe the summed responses 
using these five components of measured cardiac changes. We can 
provide each individual animal with its unique fear resistance 
index. This has enabled us not only to compare the individual 



response consistency in the same animal over the different types of 
fearful events, but also use the cardiac dynamic-based metrics for 
comparative studies across different animals. It is well known that 
there is remarkable inter-individual variability in human's stress 
responses [52]. Complex interactions between genes and devel- 
opmental experience could cause the variation in humans [52]. 
The present study has investigated the individual differences of 
cardiac responses to fearful stimuli in isogenic mice. Our result 
showed each individual mouse had consistent cardiac responses to 
different fearful stimuli. Furthermore, we noted significant 
differences among our individual wild-t^pe mice (C57BL/6J). 
Although the inbred strains of mice, such as C57BL/6J mice are 
genetically homogeneous, the behavior in mice may result not only 
from genetics but also from environmental influences [53]. Several 
studies have showed the individual differences of behavior in the 
isogenic mice [54-57]. The prenatal period and post-natal period 
experiences, such as intrauterine position, nutrition in utero, 
maternal stress, and social status may exert certain influences 
[58,59]. This suggests the possible epigenetic factors may 
contribute to such variability of cardiac responses to fearful 
stimuli in C57BL/6J mice. It is conceivable that the fear resistance 
index can be quite useful in future characterizations of neural 
responses in the mouse brain, especially in relationship with the 
processing of real-time fear memory traces associated with air puff, 
free fall and earthquake experiences [28-30] . 

In conclusion, we have shown the rapid increases in HR and 
great reduction in HRV as mice were subjected to fearful events 
such as earthquake, free fall, and air puff. By systematically 
characterizing five distinct components related to changes in HR 
and HRV, we show that the duration of the maximum plateau 
phase and HR recovery speed are quite sensitive to the habituation 
effects. Moreover, based on the measurement of the fear resistance 
index, we can compare the summed cardiac responses across 
different fearful events and across different animals. The isogenic 
mice (C57BL/6J) showed consistent differences on inter-individual 
fear response magnitudes. In addition, it is generally accepted that 
in affective disorders altered HR dynamics will increase the risk of 
cardiac mortaUty [49]. Therefore, the fear resistance index may 
open a door to measure personality traits or individual differences 
in stress-susceptibility in both wild-type mice and post-traumatic 
stress disorder (PTSD) models. 

Materials and Methods 

Ethics Statement 

AU animal work described in this study was carried out 

according to the National Institutes of Health Guide for Care 
and Use of Laboratory Animals, and was approved by the 
Institutional Animal Care and Use Committee of Georgia Regents 
University and Banna Biomedical Research Institute. 

Subjects 

A total of seventeen adult male C57BL/6J mice were used for 
ECG recording. Mice were individually housed on a 1 2: 1 2-li light- 
dark cycle (Lights on at 8:00 a.m.) and had access to food and 
water ad libitum. Data from thirteen mice in which the ECG signals 
were stably recorded were analyzed here. 

Animal Surgery and Electrocardiogram (ECG) Recording 

Electrocardiogram (ECG) measurements were employed on 
freely behaving mice, similar to that previously described [42]. 
Briefly, mice were isolated from standard cages and singly housed 
in the round home cages (48 cm in diameter, 42 cm in height) one 
week before surgery and handled twice a day (5 min each session) 
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to minimize the potential stress from the experimenter. On the day 
of surgery, the mouse was anesthetized with i.p. injection of 
Ketamine/Domitor (60/0.5 mg/kg). With the mouse placed 
supine on a heating pad, a single pair of insulated electrodes 
penetrated subcutaneously from the back of the neck to the chest 
in the lead II configuration [60] . The positive electrode was placed 
in the left abdomen below the diaphragm and anchored to the 
underlying peritoneal tissue, while the negative electrode was 
placed in the right upper chest and anchored to the pectoral 
muscle. The incisions were then carefully sutured. The mouse was 
a\v()k(;n with an i.p. injection of Antisedan (2.5 mg/kg), and then 
returned to the home cage to allow for fuUy recovery (at least five 
days) before any of the behavioral tests were introduced. 

The ECG electrodes were pre-attached to a miniature 
connector (Omnetics). The latter was then connected to Multi- 
channel Acquisition Processor system (Plexon Inc., Dallas, TX) via 
ultra-fme 36 gauge wires (0.127 mm in diameter, Plexon Inc., 
Dallas, TX). A small helium-filled mylar balloon was tied to the 
middle of the wires to balance the weight. The wires were long 
enough (230 cm) so that the animal was able to move freely within 
the home cage and experiment chambers. The ECG signals 
(filtered at 0.7-300 Hz, digitized at 5 kHz) were monitored and 
recorded by using the Plexon Sort Client software, as well as the 
mice behaviors by using the Plexon CinePlex video recording and 
tracking system. 

Behavioral Paradigms 

Three different types of fearful stimuli were designed and 
performed as previously described [28] . (1) A sudden air puff to the 
animal's back (termed air puff, 10 psi, 0.5 s). The air tube (a vinyl 
tube, I.D. 1.5 mm) was bundled with ECG recording cable. 
Therefore, the air puff was consistendy delivered to the animal's 
back, when the animal was freely behaving in the home cage. (2) 
An abrupt drop of the animal (termed free fall, 1 3 cm high) inside 
a small elevator (a square chamber, 10x10x15 cm). (3) An 
unexpected earthquake-like shaking to the animal (termed 
earthquake, 3000 rpm, 0.5 s) inside a circular chamber (12.5 cm 
in diameter, 15 cm in height) fixed on a vortex mixer. To maintain 
the consistency of stimulus inputs and yet minimize possible 
prediction of upcoming stimuli, the fearful events were triggered 
using a computer. The air valve (ASCO Valve, Inc., SC 
8202A206V), suspension system (Magnetic Sensor Systems, Series 
S-20-125), and vortex machine (Thermolyne Maxi Mix II Type 
37600 Mixer) were precisely controlled in timing and intensity by 
the computer programs. The timestamps of the fearful events were 
recorded along with the ECG signal. To examine potential 
habituation effects, we delivered each type of fearful events to mice 
for seven times at randomized intervals varying between 2 and 
5 min. The experiment sequence was as follows: earthquake, free 
fall, and air puff. The interval between sessions was 1 h. In order 
to minimize the effect of novel environment on HR [42] , mice 
were habituated in both the free fall and earthquake chambers for 
1 h per day, and three days in total. On the day of experiments, 
mice were further allowed to be habituated to the chambers for at 
least 15 min before free fall and earthquake stimuli started to 
deliver. 

In a control experiment, a second group of mice were subjected 
to neutral tone (30 s, 85 dB, 5000 Hz) in the home cage seven 
times with 2-5 min randomized intervals. 

Analysis of Cardiac Response Dynamics 

All the ECG analyses were conducted in Matlab by using 
custom-written codes. The R-wave maxima of ECG signals were 
obtained by peak detection algorithm to calculate the R-R 



intervals, and the R-R interval time series were converted into 
instant HR time series for analyses. Upon a fearful stimulus, the 
HR first increased dramatically (rising phase), and then entered 
into a maximum plateau phase during which HR remained high 
and HRV greatiy reduced. At last, the HR gradually recovered to 
basal levels (recovery phase). Five HR model parameters (rising 
speed, plateau duration, HR of maximum plateau, CV of 
maximum plateau and recovery speed) were derived from these 
three phases. 

For the maximum plateau detection, the point of maximal slope 
variation was chosen as the start of the maximum plateau. The 
momentary slope variation of the nth point (the 1st point was the 
onset of the stimulus) was calculated by the angle of two lines as 
follows: AS„= ta.n6„ = {k„2 — k„i) /{I +k„i-k„2), where k„i repre- 
sents the slope of a line fitted to data points {HR„-4, HR^^^, HR„_2, 
HR„ J, HR„}, and k„2 represents the slope of a line fitted to data 
points {HR„, HR„+i, HR^2, HR„+.i, HR^4}. The linear least 
squares fitting method was used in the line fitting. Beginning with 
the start of the maximum plateau, the accumulative standard 
deviation was calculated with 1 -point increment. The accumula- 
tive standard deviation as a function of time was derived. The 
accumulative standard deviation increased dramatically in the first 
a few points and then became steady. The sudden increase point 
after the steady state was defined as the end of the maximum 
plateau. The plateau duration was the time between the start and 
end of the maximum plateau. The HR of maximum plateau was 
the mean of instant HR during the maximum plateau phase. The 
CV of maximum plateau was based on the ratio of the standard 
deviation of R-R intervals ff to the mean of R-R intervals fi during 

the maximum plateau phase: CV= — x 100%. 

The rising phase was defined as the time between the onset of 
the stimulus and the start of the maximum plateau. The rising 
speed of instant HR was defined as the ratio of the elevation in HR 
to rising time. The elevation in HR was the difference between the 
instant HR of the first point of maximum plateau and the mean of 
instant HR during baseline period (ten seconds time period before 
the stimulus). 

The recovery phase was specified as the time between the end of 
the maximum plateau and the time point when HR recovered to 
1.1 times of the baseline. The mean of instant HR was calculated 
using a 5-s bin with a sliding window of 1 -point after the 
maximum plateau. The recovery speed of instant HR was 
calculated by the ratio of the reduction in HR to recovery time. 

We used the Poincare plot analysis to visualize the dynamics of 
HRV. Each heartbeat interval, R-R intervain was plotted on the X- 
axis against the subsequent heartbeat interval, R-R intervaln+i on 
the Y-axis. 

Root Mean Square of the Successive Differences (RMSSD) is 
one of time-domain tools used to assess HRV, the successive 
differences being neighboring R-R intervals: 



RMSSD = 



A 



1 



N-l 



{{R-R)i+i-{R-R)f], where {R 



R}i represents the ith R-R intervals, jV= number of R-R intervEil 
terms. 

Fear Resistance Index 

To obtain the summed cardiac responses from a given mouse in 
response to a given fearful event, we calculated the weight of each 
of the five cardiac parameters by using the Shannon's Entropy 
Weight Method. These five parameters are as follows: (1) the rising 
speed (of the rising phase) which characterizes the speed of the 
elevation of HR upon stimulation; (2) time duration of the 
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maximum plateau; (3) HR of the maximum plateau; (4) CV of the 
maximum plateau; and (5) recovery speed (of the recovery phase). 
Supposing there are n evaluation parameters and m sample data, 

the original data matrix X = I x, - ) is formed. The main steps 

of using entropy coefiicient method to determine the weights are 
as follows: 

(1) The standardization of original data. The raw data are 
standardized to eliminate anomalies with different measure- 
ment units, scales and positive-negative orientation. This 
process transforms different scales, units and orientation 
among various parameters into common measurable units to 
allow for comparisons of different parameters. The parameter 
playing a positive role: x,y = .\:,y /maxXy, where maxXj = 
max(x,y|; = 1,2,. ..,»?). The parameter playing a negative role: 

Xij = mmx'j^ x'-j, where minxj= min(x,y|z= 1,2,. Then 
setPij= -i^, j=l,-,n. 

i=\ 

(2) Computation of entropy. The entropy of the jth parameter 

m 

was defined as follow: Hj= —k^^Pij In Py j=\,2,---,n, 

! = 1 

where k is the entropy constant and is equal to 1 /In (m). 

(3) Evaluation of parameter weight. The degree of diversification 
is calculated by dj = 1 — H/. The weight ofj'th parameter is 

n 

Wj = djl Y^dj , which indicated the degree of importance of 

7=1 

parameter J. 

(4) Comprehensive evaluation value. The comprehensive evalu- 

n 

ation value can be derived as Uj = ^VjPij. 

,/ = l ^ 

Here, the comprehensive evaluation index, namely fear 
resistance index of different mice for three types of fearful stimuli, 
was calculated respectively. Take earthquake for example, seven 
sample data points from the first trial of seven mice and five 
parameters (rising speed, plateau duration, HR of the maximum 
plateau, CV of the maximum plateau and recovery speed) were 
used to form the evaluation matrix. The greater fear resistance 
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